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ly^ . Abstract 



A consistent combination of quantum geometry effects rules out a large class of 
models of loop quantum cosmology and their critical densities as they have been 
O^, used in the recent literature. In particular, the critical density at which an isotropic 

^ \ universe filled with a free, massless scalar field would bounce must be well below the 

Planck density. In the presence of anisotropy, no model of the Schwarzschild black 
hole interior analyzed so far is consistent. 



Q>^ ■ Aside from detailed technical constructions, one of the main achievements of loop quan- 

. tum gravity [H El E] is to provide a framework for discrete dynamical geometries. Loop 
quantum cosmology [1] takes these ingredients and applies them to expanding universes or 
black hole models. In this way, effects of quantum physics as well as quantum geometry 
can be explored in detail. This is a specific realization of a general issue which has been 
discussed recurrently in the context of discrete models of gravity O EJ [TJ [H [9] . Loop quan- 
tum gravity makes many of the older considerations more specific, but since neither its 
complete form nor the precise transition to loop quantum cosmology has been hammered 
out as of now, the dynamics of loop quantum cosmology cannot be unique. Current ig- 
I norance must be parameterized so that at least qualitative implications can be found. In 
particular, all corrections from quantum physics and quantum geometry must be analyzed 
before reliable and robust conclusions can be drawn. Here we consider quantum geome- 
try in regimes where its implications are dominant over genuine quantum corrections, and 
provide a consistent combination of its two main effects: holonomy corrections but also 
inverse volume corrections which have often been ignored. As we will see, this combination 
leads to tight consistency conditions which rule out parameter choices made so far in cases 
where only holonomy corrections were considered. Quantum back-reaction effects, which 
still remain to be fully derived following the methods of [TOl [TT], will not be required for 
our analysis. 

We start with a discussion of isotropic models, which classically have only one invariant 
scale: the Hubble distance 7i~^ = a/a (a dot meaning a derivative by proper time). If we 
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consider a fixed region V in our isotropic space, a second scale a^Vo arises which is inde- 
pendent of coordinates but depends on the coordinate size Vq of the region chosen. These 
two scales provide the classical canonical pair, {7i, Voa^} = 4:TtG with the gravitational 
constant G. 

Quantum geometry effects are expected to arise if the Hubble scale is Planckian, 
Ti.^^ ~ ip. Technically, this is realized in loop quantum gravity because holonomies, 
which are nonlinear functions of a and a, replace the classical Hubble parameter in the 
gravitational Hamiltonian [12]. But also the intrinsic spatial geometry, which depends on 
a but not on a, implies corrections when all of space is made of discrete patches whose 
size is near Planckian. (By "patch" we will mean the smallest building block of a discrete 
geometry.) This arises technically because inverse patch sizes appear in Hamiltonians, and 
their quantizations differ from classical values for small enough patches [T3l[Tl]. However, 
in contrast to the invariant scale a /a, there is no invariant measure for the spatial geometry 
in a purely isotropic setting: the scale factor a is coordinate dependent while the coordi- 
nate independent size VoO^ depends on the region V chosen. This is the reason why inverse 
volume corrections have not been fully realized in minisuperspace models so far, and most 
often were ignored. 

The patch size refers to an underlying discrete state giving rise to the expanding uni- 
verse. To include this effect in a pure minisuperspace model, extra input is thus required 
from the full theory, which must guide the construction of homogeneous models. Fortu- 
nately, this is possible in a heuristic mean-field picture [15] where a mean-field describing 
the underlying discrete structure is used. As always with mean-field approximations, the 
precise way in which the field arises from a fundamental theory may be very difficult to de- 
rive. In lieu of a detailed derivation, however, there are often consistency conditions which 
can help to reduce the a-priori freedom involved in the choice of mean-field behaviors. This 
turns out to be the case for loop quantum cosmology, too. 

We then have a patchwork structure describing the isotropic classical geometry in an 
atomic manner. Instead of the classical size a^Vo, we can write J\fv where v is the mean 
size of a patch and A/" is their number in the region V. Since both products are supposed 
to give the size of the same region, they must equal each other: 

Xit)v{t)=a{tfVo. (1) 

In this equation, we have already accounted for the main result from loop quantum gravity 
employed here: the fact that any dynamics proposed for the full theory [THl [131 [IT] must 
change the values of discrete contributions to volume; v(t) must depend on time as nodes 
representing the atomic volume elements are being excited. In general, also the number of 
patches J^(t) is a function of time, and the product of both functions provides the time 
dependence of the scale factor a(t). (For the current purposes it is irrelevant what time 
variable is used. In explicit constructions, this could be an internal time such as a matter 
field.) 

As an important consequence, we see that the expansion of an isotropic universe in 
loop quantum cosmology is described by two time-dependent functions, Af(t) and v{t), 
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rather than one classical function a(t). Thus, there is more freedom in the underlying 
dynamics which is related to the refinement behavior of an underlying discrete state. By 
a pure minisuperspace quantization, which starts from the classical a and turns it into an 
operator, this freedom cannot be constrained but rather emerges in the form of quantization 
ambiguities in the reduced Hamiltonian. The detailed form of the dynamical change of 
A/" and V can in principle be derived from the full theory, or be restricted by analyzing 
the implied phenomenology of models. (It is not uncommon that additional parameters 
and functions not seen classically arise from a quantum or microscopic treatment. In 
cosmological models, for instance, a similar behavior has been observed in [18j.) 

Quantum geometry is discrete but local, and thus corrections to the classical geometry 
and its dynamics depend only on the microscopic patch size v rather than the macroscopic 
number A/" of patches. (From a classical perspective, this behavior may not seem local if 
a discrete patch such as the geometry given by a single spin network vertex is viewed as 
a representation of the continuum distribution in a neighborhood. But quantum geom- 
etry is local in the sense that geometrical operators only refer to sums over single patch 
contributions, a behavior which we will call "patch-local" in what follows.) This already 
implies that all consistent implementations of quantum corrections must be independent 
of the size Vq of the region, as this is a parameter which determines A/" but does not affect 

V. 

We first address holonomy corrections which are easier to construct, although also 
here extra input is needed compared to the original and direct quantization used in loop 
quantum cosmology [19]. Holonomies are non-linear functions of a connection which in 
isotropic models is simply proportional to d: c = 7(2 with the Barbero-Immirzi parameter 
7 I20I [21]. In a patch-local holonomy, the connection appears in a line integral along a 
curve of the (linear) patch size Iq. In an isotropic setting, holonomies are thus functions of 
-^oc = '~fLa/a where L = i^a is the geometrical and coordinate independent length of the 
curve (or its co- moving size). 

The number of patches of coordinate size £9 in a region of size Vq is given by A/" = Vq/^q- 
Together with ([1]), this implies that L = ioa = {Vo/J\fy^^a = v^^^ is directly related to 
the microscopic patch size v. Thus, holonomy corrections, which arise from deviations of 
non-linear functions such as sin(7Ld/a) from the linear 'jLa/a, are independent of the size 
Vq, as required. Below we will see that this is also true for inverse volume corrections. 

Before discussing inverse triad corrections and their mutual consistency with holonomy 
corrections, we emphasize the dynamical role of L{t) = v{tY^^. In the original formulation 
of loop quantum cosmology [121 [22], it was assumed that £0 is constant as the universe 
expands. This is the most straightforward assumption in a pure minisuperspace construc- 
tion, but it was clear that this cannot describe long phases of expansion: a constant 
£0 = {Vo/J^y^^ implies that the number of patches in the co-moving region V is constant. 
As the geometrical size of this region is expanded by the scale factor, a constant number 
of patches means that patch sizes are expanded to macroscopic values. This would clearly 
be in conflict with a large universe free of noticeable quantum geometry effects. (The cor- 
respondence between constant io and a constant number of patches is also demonstrated 
by the two-patch model of |23j.) 
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To address this problem, an improvised version of the dynamics was proposed in [23], 
based on the ad-hoc assumption that low-lying eigenvalues of the full area spectrum should 
determine quantization parameters of minisuperspace models. This assumption is ad-hoc 
because the area operator plays no role in full constructions of Hamiltonian constraints 
or their solutions. One could think that it somehow arises in a gauge-fixing procedure 
which leads one from the full theory to an isotropic model, but no such construction is 
known. (The use of the area operator, which is not a physical observable, was criticized 
in [25]. However, the area is used when constructing the constraint and before solving it; 
non-observables can certainly appear in this step just as a/a, which is not a constant of 
motion, appears in the classical constraint.) Despite of these shortcomings, the procedure 
implied that it is not Eq which is constant but L = i^a = v^^^. This means that the 
microscopic patch size is constant in this version, and correspondingly the number of 
patches increases in an expanding universe. In this way, strong quantum geometry effects 
in a large universe can be avoided for suitable parameter choices. However, the values 
proposed in [23], motivated by the smallest non-zero area eigenvalues, turn out to be in 
conflict with inverse volume corrections, to which we turn now. 

Inverse volume corrections arise because the patch volume v, when quantized, becomes 
an operator with discrete spectrum containing zero [261 [2Z]- Such an operator does not 
have a densely defined inverse, but an inverse patch volume is needed in most Hamiltonians. 
One can construct operators with the correct classical limit given by the inverse based on 
identities such as 

for the canonical pair c = Vq''^c = VQ^^'yd and \p\ = Vq^^o^, {c,p} = Sir^yG/S, under- 
lying isotropic loop quantum cosmology. For < r < 2/3, the right hand side gives an 
inverse power of the scale factor, while no inverse is needed on the left. This gives rise to 
well-defined quantizations, which are part of a general procedure [13], and whose specific 
constructions in isotropic models all follow the derivations in [28] to which we refer for 
further details. 

Since the quantized holonomy exp{iiQc/VQ^^) becomes a shift operator in p and the 
Poisson bracket in ([2]) becomes a commutator, the quantization of inverse volume pro- 
vides a discrete approximation of the derivative involved in the classical Poisson bracket. 
The derivative d|pp'"/^/dp on the right hand side of ([2]) is then replaced by a difference 
{\p + Apl^*^/^ — \p ~ Ap|^''"/^)/2Ap where Ap depends on the precise implementation of the 
quantization scheme and is typically related to the Planck length. (Numerical studies of 
the volume spectrum in the full theory, whose results will shed light on which precise values 
are to be expected, are on-going [291 EQl [311 [32].) Deviations of the difference from the 
derivative give rise to inverse volume corrections. 

The form of Ap is important to demonstrate the consistency of inverse volume correc- 
tions, i.e. independence of Vq, and to estimate their size. Inverse volume corrections as 
they have been used in most isotropic models so far correspond to using the total volume 
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V = a^Vo of the region V and approximating its derivative by differences of the form V±ip. 
This has two consequences: (i) inverse volume corrections of this form are extremely tiny 
when the total volume is much larger than the Planck volume, and (ii) such corrections 
depend on the region V and its coordinate size Vq because V does. The first property 
led [21] to conclude that inverse volume corrections can safely be ignored compared to 
holonomy corrections, but the Vo-dependence remained unexplained; in fact, it points to 
an inconsistency in the derivations used which cannot simply be eliminated by ignoring 
inverse volume corrections. 

In the present context one can see a third problem with such a treatment of inverse 
volume corrections: if quantizations of the inverse volume depend on the total volume of a 
region, resulting operators cannot be patch-local. Inverse volume operators appearing e.g. 
in energy densities would depend on the full volume of a region V which can be macroscopic. 
Although this may not appear as a problem in an isotropic context where homogeneity 
identifies local and global properties, it is clear that such quantizations cannot arise in a 
model consistently related to what we know from the full theory. As already mentioned, full 
Hamiltonians have contributions which are patch-local and depend on the patch-volume 

V but not on the number of patches A/" or the total volume V. This has two important 
consequences: (i) consistent implementations of inverse volume corrections do not depend 
on Vo since v does not, and (ii) they are much larger than naively expected because 

V is much smaller than V and differences v ± ip deviate much more from infinitesimal 
displacements of v than V ± ip does from displacements dV^ of the large V. That inverse 
volume corrections are artificially suppressed in direct minisuperspace quantizations has 
also been discussed in the appendix of [33]. 

Consistent inverse volume corrections imply correction functions a{L/ip) at all places 
where inverse densitized triad components are used in a classical Hamiltonian, such as 
in the kinetic term of a matter Hamiltonian [Tl|. The precise functional form of a is 
subject to quantization ambiguities [Ml |35l EE] , but is restricted by anomaly cancellation 
conditions |37j. In all cases, a deviates strongly from the classical value a = 1 when 
its argument is of the order one or smaller; for large values of L/£p the classical limit is 
approached asymptotically. Importantly, correction functions a{L/ip) are to be evaluated 
at L/£p rather than the much larger and Vo-dependent V^^^/ip. There is thus no Vo- 
dependence, and quantum geometry effects from inverse volume operators are much more 
pronounced than expected. In this way, valuable consistency conditions can arise. In 
particular, for holonomy corrections to be small in classical regimes we require ioc = 
'jLd/a <C 1, while small inverse volume corrections require v^^^/ip = L/ip ^ 1. For a 
given classical geometry, we thus have an upper as well as a lower bound on L, which leaves 
only a finite range of allowed choices. 

To make this quantitative, we derive the main effect of holonomy corrections in the 
current setting, and compare parameter values with what is consistent with inverse volume 
corrections. The Hamiltonian constraint is the main place where corrections enter, and 
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holonomy corrections imply a form 



C = -^^^V^si^'(^oc/l^o'^') + \P?"P = (3) 

with a matter density p. These terms strictly refer to operators as they appear in a 
quantization of the Hamiltonian constraint, but they can also be considered as the tree- level 
approximation to the quantum constraint operator. Our considerations will not require 
deep quantum regimes as they are, e.g., encountered around classical singularities, where 
such tree-level equations would be unreliable due to the presence of additional quantum 
corrections [38[ |39] . The tree- level constraint then implies equations of motion with tree- 
level corrections from quantum geometry. We will only require the equation 



T/l/3 

p = {p, C} = sm{ioc/V,'^') cosMV^'''^ 



1^0 



which can be written as 



2 / • \ 2 

a\ I p 



a) \2p) ;^-=(^oe/vV»)(l-.„=(V/V^r)) (4) 

5^, U _ (5) 



where we used the Hamiltonian constraint ([3]) to express sin^(£oc/V"Q^''^) in terms of p. 
Defining the critical density 

we thus have the tree-level Friedmann equation [ID] 

'"^pfl-^l (7) 



a J <3 V Pcrit 

whose behavior is determined by the patch size L (see also [H] for the derivation with 
general L). The meaning of the critical density pcrit can be seen in a model where the only 
matter ingredient is a free, massless scalar: In this case the tree-level equation is an exact 
effective equation [12] and can be used even in deep quantum regimes. When p = p^rit, 
the scale factor reaches a minimum and the universe bounces. For other matter there are 
additional quantum corrections [39], and even the tree- level form can change if higher order 
holonomy corrections are included which are possible due to quantization ambiguities [13] . 
The bounce has reliably been demonstrated only for a free massless scalar, but the critical 
density can nonetheless be used more generally as a measure for the size of holonomy 
corrections. 

It is hard to derive L from an inhomogeneous state in full loop quantum gravity, but 
one can easily check the mutual consistency of any value for pcrit in holonomy corrections 
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with the corresponding imphcations of L in inverse volume corrections. In the case of 
[23], which is the main specific model analyzed recently, we have a constant L which, 
by reference to the full area spectrum, gives rise to a critical density near the Planck 
density: pcrit ~ Pp = l/G£p. With ([6]) we obtain 7L ~ a/S/Stt^p and thus L/£p ~ 
■\/3/87r7~-'-. For this value we have L/£p > 1, thanks to the smallness of 7 ~ 0.24 as it 
arises from black hole entropy calculations [Ml . This is sufficient to avoid extremely 
strong inverse volume corrections. However, L/ip ~ 1.4 is not much larger than one 
and inverse volume corrections for this value are still significant. In the specific case of 
constant L as in [23] , the consequences may not be too severe because a constant L implies 
that inverse volume corrections merely amount to a constant factor in all terms of the 
constraint where an inverse volume appears. This changes some coefficients, but would 
leave the qualitative behavior untouched. However, such a parameter choice would be 
unstable: a slight deviation from a constant patch size L, as it can easily occur in general 
states, would mean that inverse volume corrections become a-dependent. Now, it would 
not only be coefficients but even the qualitative dynamical behavior which is subject to 
strong changes since, for instance, corrected energy densities p containing a{L/ip) would 
depend differently on the scale factor than they do classically. 

In a strict sense with the usual value of the critical density, we have disproved the 
improved dynamics of [23]. This dynamics can be consistent with weak inverse volume 
corrections only if the critical density is much lower than Planckian, and the patch size 
much larger. For L larger than ip by one order of magnitude, for instance, the critical 
density could at most be about 2% of the Planck density. 

A constant patch size f or a constant L is only a special case, for which the number 
of patches is proportional to the volume, oc a^. In general, however, both v and A/" 
can be functions of a if the scale factor is used as internal time to measure the change 
of discrete geometry. As usually, it is often convenient to assume power-law behaviors at 
least for certain phases of the universe evolution, which we parameterize by a power x in 
A/" = Aoct"^^. The patch size then behaves as 

V = a^Vo/AT = 0=^(^+2" Vo/ATo . (8) 

Notice that the parameter A/q in general depends on the size of the region V as well as 
coordinates in order to make Af coordinate independent and proportional to Vq; it is not 
simply a numerical constant. 

In this parameterization, we have x < if the number of patches A^ is increasing, 
and X > —1/2 if the patch size v is increasing. Since there is a lower limit for both Af 
and V, values x > or x < — 1/2 cannot be realized forever. Thus, if a uniform power- 
law model is used, as it is for instance done in solvable models [^ |36[ 133 I3H] for a free 
massless scalar, one has to restrict the power to the generic range —1/2 < x < 0. The 
limiting values are possible, but would not be generic as they would require either no 
change in the patch size to occur [x = —1/2) or a constant number of patches [x = 0). 
The former choice, which gives a constant L and thus the dynamics of [24], is impossible 
in full constructions of loop quantum gravity while the latter is not consistent with large- 
scale semiclassicality. In general, power laws can only be assumed for certain periods of 



7 



time, as it is commonly done for energy densities in different phases of the universe, and 
occasionally a power x may even fall outside the range —1/2 < x < 0. Different power- 
laws can then follow each other in much the same way as, e.g., a dust-dominated phase 
with one power-law for energy density would follow a radiation-dominated phase with a 
different power- law. Complete functions J\f{a) and v{a) incorporating these different phases 
could only be derived from a full solution of an inhomogeneous state which is difficult, but 
phenomenological constructions are already possible. See e.g. [HI |50l EH |52] for early- 
universe restrictions on x and L. As control over the theory increases, one can expect tight 
consistency conditions to arise in this way. 

To restrict general choices even more strongly, anisotropic models are valuable. Effects 
of lattice refinements, as first analyzed in this context in [33], now allow more freedom 
because parameters in general depend not only on the volume but also on extensions 
in different directions. The interior of a Schwarzschild black hole, for instance, is of a 
Kantowski-Sachs geometry with two independent scales to determine the spatial geometry: 

ds2 = l^da;2 + |p2|dl]2. (9) 

There are two densitized triad components pi and p2, which provide the volume V = 
AitLq\pi\^,/\p^ where Lq is a coordinate length parameter to select a finite region whose 
angular extension has been fixed by using the full sphere size 47r; for details see [53] or 
[5l] . When triad components and thus the volume vanish, we have either the black hole 
singularity {p2 = 0) or the horizon (pi = 0). 

In this context, we can write a more general equation to relate microscopic quantities 
to the classical geometry: 



Af,{t)L{t)Af2{t)A{t) = 47rLobi(t)|/hM (10) 

where A/i is the number of patches in the x-direction, with patch extension L in this 
direction, and N'2 is the number of patches in the spheres spanned by the angular directions, 
with patch area A in these directions. (Spherical symmetry implies that we do not have to 
distinguish the two linear extensions in angular directions, which we thus combine in the 
angular area.) 

The volume relation fllOl) then factorizes in two identities 

M(t)L(t) = Lo^ii|L , ^^mit) = A7r\p2{t)\ (11) 

using the sizes following from the spatial Kantowski-Sachs metric in triad variables. 
To specify the refinement fully, we have to determine two functions such as J^i{pi,P2) 
and J^2{Pi,P2)'i just the volume dependence of the total number A/1A/2 of patches is not 
sufficient. A particular case, which provides constant patch sizes and shapes, is obtained 
when L and A are constant and thus A/i oc \pi\/ a/|p2| and A/2 oc \p2\. This case was 
introduced in [53] and shown to have stable and near classical behavior in a large region of 
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the phase space. It is a special case of models where the total patch number is proportional 
to volume: MiN'2 oc \pi \ ^J\p2\ oc V. 

More generally, one can assume a volume dependence MiN'2 oc which still leaves 
different choices for the anisotropic behavior of the patches. What is important in the 
present context is that none of these power-law cases can be fully consistent in classical 
regimes: The total volume vanishes not only at the singularity, where strong discreteness 
effects would not be unexpected, but also at the horizon where pi = 0. For a consistent 
model, which behaves semiclassically near the horizon as it should be the case for massive 
black holes, one has to use more general functions than power laws to parameterize the 
refinement behavior. In particular, A/1A/2 must not vanish when pi = 0. Despite many at- 
tempts, no fully consistent effective Schwarzschild interior space-time has been constructed 
yet, even if one restricts oneself to regions far from classical singularities where tree- level 
equations are reliable and quantum back-reaction effects can be ignored. Anisotropic mod- 
els thus provide important test cases which are restrictive enough to constrain lattice 
refinement models further. 

To conclude, we have seen that a consistent treatment of loop quantum cosmology, with 
a robust connection to constructions in the full theory, requires input from the behavior 
of an underlying microscopic state. If this is done correctly, the balance of quantum 
corrections changes such that several choices made so far can be ruled out because they 
would imply strong quantum geometry corrections. (In a pure minisuperspace setting, 
models with suppressed quantum geometry corrections can be constructed consistently and 
rather straightforwardly. But they do not consistently model the full theory.) In particular, 
critical densities must be sufficiently smaller than Planckian, which enlarges holonomy 
corrections; this may be of interest for potential observations in cosmology because effects, 
such as those in the tensor mode spectrum [^Sl ESI EZl EHl ES], would be enlarged. This 
would have no effect on existing analyses of inverse volume corrections for perturbative 
inhomogeneities [601 133 EI] where consistent versions were already used. Since it has 
become quite fashionable to boldly extend tree-level equations such as ([7]) and its analogs 
in black hole models well beyond their proven range of validity and deeply into putative 
quantum regimes around classical singularities, one should note that this is not done here. 
To rule out a Planckian critical density we do not have to consider the dynamics near a 
classical singularity, nor do we have to consider long stretches of evolution which would 
both give rise to quantum back-reaction in non-solvable models. For a Planckian critical 
density, consistent inverse volume corrections would be significant at all times even if 
matter densities are small. In this way, reliable upper bounds on the critical density arise. 
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